Spatial Analysis

Load data

Maps data

We use three different shapefiles for the continental U.S. land mass, the State waters of maine, new hampshire, massachusets, connecticut, rhode island, new york, new jersey, delaware, maryland, virginia, north carolina and the North East EEZ.

1.- Land shapefile; covers the US land territory for visualization. Data provided from the map package.

2.- State waters; covers the state waters of the NE US states. Data from data.gov.

License: No license information was provided. If this work was prepared by an officer or employee of the United States government as part of that person’s official duties it is considered a U.S. Government Work.

3.- EEZ shapefile; Used the Sea Around Us shapefle updated to June 2016.

unique(land_sf$ID)

Spatial component

Create base shapefile for computation

As a first step we need to divide the NE US EEZ among the different states. For that we expanded state waters up to the 200 nautical miles to then estimate the percentage that each expanded-state-waters occupied. Note that in all cases these areas overlapped and percentages accounted for it. We did this by following these steps:

• 1. Expand state waters using a buffer
• 2. Grid that buffer on a 0.3 resolution
• 3. Crop the buffer to the EEZ

1. Expand Spatial regions (a.k.a buffers)

We set a buffer of 410000 m (410 km, ~ 221 nm) that overshoots the EEZ a bit, but is eventually cropped

# Buffer state waters
state_bf = st_buffer(state_sf, 410000) %>% 
  st_transform(4326) %>% # to match shape
  st_set_crs(4326)

# ------------------------------ #
# Testing and visualizing the buffer 
# ------------------------------ #

state_bf_plot <- ggplot() +
  geom_sf(data = eez_sf, aes(), fill = NA) +
  geom_sf(data = state_bf, aes(color = state), fill = NA)
#   

# ggsave("../Results/Partial/state_waters_buffer.png",state_bf_plot)

# ------------------------------ #

state_bf_plot

2. Expand grid within buffer

Here we expand a grid within the buffer so we can estimate the proportion of each state

2.1 Merge grid and buffers

Once we have a gridded area, we converted the grid to a sf so we can merge it with the buffered states and finally filter out everything outside the states polygon

# ---------------- #
# Convert grid to sf
# ---------------- #
grid_sf <- st_as_sf(ne_grid,
             coords = c("lon", "lat"),
             crs = 4326) %>% 
  st_join(state_bf) %>% 
  filter(!is.na(state))

# Create data frame for future computations
# Note, will be used in next chunk
grid_bf_dt <- as.data.frame(grid_sf) %>%
    select(index,state)
    # group_by(state) %>% 
    # summarise(length(index))


# ---------------- #
# [TEST] 
# Visualization of grid
# ---------------- #

gridExtra::grid.arrange(
  # Overall (overlapping) position
  ggplot(grid_sf) +
    geom_sf(aes(color = state), alpha = 0.3) +
    geom_sf(data = eez_sf, aes(),fill =NA) + 
    theme(legend.position = ""),
  # Showing each state separatley
  ggplot() +
    geom_sf(data = grid_sf, aes(color = state),size = 0.1, alpha = 0.3) +
    facet_wrap(~state) +
    theme(legend.position = "top"),
  nrow = 1)

3. Crop buffers to EEZ

Finally, we crop the grided buffers to within the EEZ to capture the actual water space.

Note: This step takes quite a while because of the size of the EEZ shapefile. No, you cannot use st_simplify() here

Interpolation rutine

Here we interpolate the survey data within the previously created grid.

• We removed cases where wtcpue < 0

Function for interpolation

This is the main function used to interpolate the data per year. It follows a Triangular Irregular Surface method using the interp::interp() function.

Ocean adapt data

• Using only the Northeast US Fall and Spring bottom trawl survey data for now

Control Pannel

Load required data. Note that some of it has been previously created

# Load grid df
grid_eez_df <- my_path("D","Spatial","grid_eez_df.csv", read = T)

# Run interpolation for all years
years <- seq(1973,2019,1)

# taxa
spp <- "Centropristis striata"

# regions
regions <- c("Northeast US Fall" , "Northeast US Spring")

Run routine

Extrapolated data

index	state	year	region	spp	lon	lat	value
73	north carolina	1973	Northeast US Fall	Centropristis striata	-79.85582	30.94383	NA
74	north carolina	1973	Northeast US Fall	Centropristis striata	-79.55582	30.94383	NA
75	north carolina	1973	Northeast US Fall	Centropristis striata	-79.25582	30.94383	NA
76	north carolina	1973	Northeast US Fall	Centropristis striata	-78.95582	30.94383	NA
77	north carolina	1973	Northeast US Fall	Centropristis striata	-78.65582	30.94383	NA
78	north carolina	1973	Northeast US Fall	Centropristis striata	-78.35582	30.94383	NA
79	north carolina	1973	Northeast US Fall	Centropristis striata	-78.05582	30.94383	NA
80	north carolina	1973	Northeast US Fall	Centropristis striata	-77.75582	30.94383	NA
81	north carolina	1973	Northeast US Fall	Centropristis striata	-77.45582	30.94383	NA
82	north carolina	1973	Northeast US Fall	Centropristis striata	-77.15582	30.94383	NA

Results

Points

Map

Area plot

This graph shows the proportion of each State over time

Area plot (running)

This graph shows the proportion of each State smoothed over a 10 years running mean. It helps seeing the trends better (I Think…)